Thermotropic Biaxial Nematics: Spontaneous or Field-Stabilized? 
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An intermediate nematic phase is proposed for the interpretation of recent experimental results 
on phase biaxiality in bent-core nematics. The phase is macroscopically uniaxial but consists of 
microscopic biaxial, and possibly polar, domains. On applying an electric field the phase exhibits 
substantial macroscopic biaxial ordering resulting from the collective alignment of the domains. A 
phenomenological theory is developed for the molecular order in this phase and for its transitions 
to purely uniaxial and to spontaneously biaxial nematic phases. 

PACS numbers: 61.30.Cz, 61.30.Gd 



Since their theoretical prediction, nearly four decades 
ago fH, biaxial nematics have been a constant challenge 
in liquid crystal (LC) research 0, S H They differ 
from the common, uniaxial, nematics in that they ex- 
hibit additional orientational order of the molecules along 
a second macroscopic direction, the "biaxial" or "short" 
axis m, perpendicular to the primary nematic director 
n. The expectation that the response of m to an ap- 
plied electric field could be much faster than that of n 
has been sustaining a constant practical interest in low 
molecular mass biaxial thermotropic nematic LCs (i.e. 
the biaxial analogues of the conventional nematic LCs 
used in electro-optic applications). However, it was only 
recently that strong experimental evidence has been pro- 
duced for the discovery of such biaxial nematics, first in 
bent-core systems 0, 0] and shortly afterwards in later- 
ally substituted tetrapode nematogens.jl] 

Subsequent electro-optic switching experiments Q on 
the bent-core biaxial nematics demonstrated that the re- 
sponse of the m axis to an applied field is indeed much 
faster than that of the n director. Interestingly, the inter- 
pretation of these switching experiments suggests (a) the 
existence of a high temperature uniaxial nematic phase 
with practically no biaxial response to an applied electric 
field and (b) a transition to a low temperature nematic 
phase which is optically uniaxial and can be brought to 
a biaxial state by applying an electric field perpendicu- 
lar to n. As the electric field strengths involved (a few 
V/^m) are clearly too low to produce a substantial effect 
directly on the orientations of individual molecules, the 
field-induced biaxial state is attributed to the preferential 
alignment of the m axes of pre-existing biaxial molecular 
aggregates (domains or clusters) which, in the absence of 
an applied field, are randomly distributed about n. Elec- 
tric fields of similar strength have been used to switch the 
m director in the XRD experiments of ref Q • 

The spinning-sample NMR experiments in refi do not 
involve electric fields but there, the strong external mag- 
netic field, when not coUinear with the n director, would 
orient the m axis (here identified with the direction of 
smallest diamagnetic susceptibility of the phase) perpen- 
dicular to the plane formed by n and the field direc- 



tion. Consequently, it is possible that the static sample 
consists of biaxial domains with their m axes randomly 
distributed about a common n director (parallel to the 
magnetic field) and that a macroscopic alignment of the 
domain m axes results from spinning the sample about 
an axis perpendicular to the magnetic field. 

Recent atomistic simulations of nematics made of bent- 
core molecules indicate the existence of biaxial do- 
mains. In general, such domains exhibit dielectric as well 
as diamagnetic biaxiality and would therefore be readily 
oriented by external fields. In fact, the domains found 
in ref Q exhibit local ferroelectric (i.e. biaxial and po- 
lar) nematic ordering. This endows the domains with 
a net electric polarisation which could also orient them 
macroscopically in an external field. Dielectric fluctua- 
tion studies by dynamic light scattering in the uniaxial 
nematic phase of bent-core liquid crystals suggest the 
formation of cybotactic (smectic-like) clusters [8| and, in 
one instance [9| , randomly oriented nanodomains of anti- 
clinic ferroelectric smectic ordering are proposed for the 
structure of an optically isotropic phase obtained on cool- 
ing from the nematic phase. Lastly, persistent cybotac- 
tic biaxial clusters are clearly identified in XRD studies 
of calamitic multipode nematics for which the low tem- 
perature phases are columnar [l^l- These considerations 
suggest that the existence of uniaxial phases both, with 
and without biaxial clusters, and the possibility of field- 
induced alignment of the latter, might be of key relevance 
to the understanding of phase biaxiality in thermotropic 
nematics. 

An alternative way to view the aligning effect of the ex- 
ternal field is in terms of the orientational fluctuations of 
m. The hypothesis that these fluctuations are extensive 
enough to destroy the spontaneous long-range alignment 
of the m axis has been often used as a possible expla- 
nation of why, in spite of the predictions from molecular 
theory, thermotropic biaxial nematics are not commonly 
observed in experiments. On the other hand, in analogy 
with what is known from the elastic continuum theory of 
uniaxial nematics 11| , the application of an external field 
would quench the low wave-vector orientational fluctua- 
tion modes of the m axis. Thus, a possible interpretation 
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of the observed field-induced transition to a biaxial state 
is through the quenching of the low-energy orientational 
fluctuation modes. However, estimates of cluster sizes 
and time scales of their reorientational motions indicate 
that a continuum treatment of the biaxial fluctuations 
may not be applicable. 

In any case, the conventional static formulation of the 
nematic phase free energy solely in terms of long range 
orientational order parameters refers to a single- domain, 
uniformly ordered system and can therefore convey nei- 
ther the cluster picture nor the continuum fluctuation 
picture of the field-induced transition to the biaxial state. 
In this letter we introduce a phenomenological descrip- 
tion that allows explicitly for non-uniformity of the biax- 
ial orientational order in a thermotropic nematic phase. 
The formulation is based on the biaxial cluster picture, 
allowing for a full range of cluster sizes, from single 
molecule to macroscopic aggregates, thus avoiding the 
inherent size-limitations of a continuum treatment. 

To identify the relevant order parameters in a phase 
with local biaxial nematic order we consider a nematic 
sample of N molecules in a volume V at temperatiire T 
with the director n perfectly aligned along the Z macro- 
scopic axis. Practically, the sample may be pictured as 
filling the space between parallel plates, with the Z axis 
defining the rubbing direction that aligns the director n 
on the plate surfaces. The X macroscopic axis is taken to 
lie on the plane of the plates and the Y to be perpendic- 
ular to that plane. For simplicity we further assume that 
the molecules are themselves perfectly aligned with their 
major axis z parallel to the primary director n. This re- 
stricts the local biaxial m axis as well as the molecular 
axes x,y on the plane defined by the macroscopic X,Y 
axes. We also assume that the sample-confining surfaces 
have no aligning influence on the m axis. We then assign 
to each molecule a traceless second rank molecular ten- 
sor iriab which in the principal molecular frame of axes 
{a, 6} = x,y,z may be taken to have the components 
—TTixx = m ~ niyy and rrizz = 0. Next, assuming that 
the sample is divided into a number R (< N) of clusters, 
each labeled by an index r and containing rir molecules, 
we may define for each such cluster a tensor 
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where A^, Br denote the principal axes in the cluster r 
and the index i runs over all the molecules in that clus- 
ter. The cluster tensor will thus have two non-vanishing 
principal components; these can be expressed in terms of 
a single quantity M^'^) = -M^j-^^ = m'^^y^. Obviously 
M'^''^ depends both on the size of the cluster, through 
the number of the molecules it comprises, and on the 
degree of biaxial ordering of these molecules. The short 
YWr axis of the cluster is taken to coincide with one of he 
principal axes, say Y^-. 



The rotational invariants associated with the molecular 
and the cluster tensors are respectively niabiTiab = 2m^ 
and M^J^g^M^J^g^ = 2{N&^f (summation over re- 
peated tensor indices is implied). By summing the in- 
dividual invariants of all the clusters one may define the 
following invariant quantity for each possible partitioning 
{r} of the sample into clusters: 



a{r} = 
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This quantity strictly vanishes if each cluster contains a 
single molecule, it reduces to a{r} — —1/{N'^ — 1) if 
all the clusters have vanishing M^''^ (i.e. if they are uni- 
axial) and it takes the highest possible value cr{r} = 1 if 
the sample consists of a single cluster containing all the N 
mesogens with their molecular axes x, y perfectly aligned 
along the macroscopic directions X, Y respectively. De- 
noting by cr{f} the largest value that can be obtained 
for any of the different possible ways of partitioning the 
sample into clusters in a given microstate, one may de- 
fine the cluster order parameter a as the ensemble aver- 
age a =< a{f} >. This parameter varies in the range 
< (T < 1 and gives the extent of biaxial ordering within 
the clusters but does not give a direct measure of long- 
range biaxiality. To describe the latter we use the follow- 
ing macroscopic, second rank and traceless, tensor in its 
principal axis frame A,B = X, Y, Z, 
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with principal values —qxx = Qyy = q, and < 1. The 

quantity q measures the extent of phase biaxiality. In the 
absence of an external filed, this biaxiality is understood 
to originate from the spontaneous collective alignment of 
the clusters short axes m,.. 

Due to the assumed restriction of the molecular z- 
axis along the macroscopic Z direction, the invariants 
that can be formed from the biaxiality tensor qAB are 
even powers of q (i.e. qABqAB = 2?^, qABqBcqcA = 
etc). The electrostatic interaction of the biaxial medium 
with an applied electric field, taken to have components 
Ey = E, Ex — Ez = 0, is conveyed, to lowest order in 
9 by a term hE^EBqAB = hE^q, where the scalar factor 
h reflects the magnitude of the molecular polarisability 
anisotropy in the x, y molecular plane. Accordingly, the 
leading terms in a phenomenological Landau-deGennes 
expansion of the free energy in terms of the order pa- 
rameters a and q will be 



F = aa+-a^ + la^ + -q^ + -q^ 



euq'^-hE'^q (4) 



Considering a as the primary order parameter, the co- 
efficient a is taken to be an increasing function of the 
temperature, exhibiting a rapid variation in the vicinity 
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of a characteristic temperature Tq at which it changes 
sign. The other coefficients 6, c, d, e, h, g are ah assumed 
to be insensitive to variations of temperature and pos- 
itive, with the exception of 6, for which both signs are 
considered. The c and d terms correspond to the entropic 
drop caused by the ordering and the e term is the ener- 
getic contribution associated with the couphng between 
the local ordering of the individual clusters (cr) and their 
collective ordering (q). 

With no applied field = 0), the free energy in eq (j4|) 
describes three possible nematic phases: (i) a "proper" 
nematic phase (iV„), in which a — q = 0; (ii) a macro- 
scopically uniaxial nematic phase ( A^l^'^^ ) formed by ran- 
domly oriented biaxial clusters, in which ct > and q — 0; 
(iii) a macroscopically biaxial phase (A^f '"^■') formed by or- 
dered biaxial clusters, in which a > and q 0. Repre- 
sentative order parameter profiles for these phases, and 
the possible phase transition sequences are depicted in 
Fig. dl]). 

For 6 < 0, the transition from the to the iV^*"^^ 
phase is of first order, Fig. (la), and is obtained when 
a reaches the value a*{= 3b^/{16g)), at which point the 
cluster order parameter undergoes a jump from ct = 
to CT = (J*{— 3|&|/(4(7)). A further transition from the 
Nu'"^^ to the Nj^'"^\ which is of second order, is obtained 
in this case as ct increases beyond a critical value CTc(= 
c/(2e)), provided that the ratio A = CTc/ct* is A > 1. The 
respective value of a at this transition is Oc ~ a*A(4 — 3A). 
If A < 1, the A^u*"^^ phase is removed from the sequence. 
Fig. (lb), and a direct, first order, Nu to Nj^'"'"' phase 
transition is obtained at a = ad{< a*), with both ct and 
q rising abruptly form to finite values ad and qd- 

For 6 > the transition from the Nu to the Nu"^^ 

and is 



phase is of second order. Fig. (lc,d), at a 
followed by a transition from the Nu"^ to the Nj^'"'^ on 
lowering a to the value a'' = — a*A(4 + 3A). The order 
of this transition is controlled by the parameters A and 
m(= 2e^/ (l&l d)), with m < 1 -I- 3A/2 defining the range of 
the second order phase transition. 

In the presence of an electric field {E ^ 0) the uniaxial 
phases and Nu'^^ acquire field-induced biaxiality, to 
which we now focus our attention, particularly for the 
case 6 < and A > 1 which is directly relevant to the 
biaxial electro-optic response and the nematic- nematic 
phase transitions observed experimentally in bent-core 
nematicsjl]. 

The dependence of the order parameters ct and q on the 
temperature function a is shown in Fig. 2 for different 
magnitudes of the applied electric field. It is apparent 
from the plots of the biaxiality order parameter g(a, K) 
that the effect of the electric field is much stronger in the 
Nu'^^ phase, where the biaxial ordering is produced by 
the alignment of the biaxial clusters, compared to the Nu 
phase where the field influences directly the orientations 
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FIG. 1: Plots of the temperature dependence of the calcu- 
lated order parameters for four representative combinations 
of the values of the expansion coefficients in eq (j4| : (a) 6 < 0, 
A = 1.2, u = 0.5. (b) & < 0, A = 0.7, u = 0.5. (c) b >0, 
A = 0.2, u = 1.5. (d) as for (c) with u = 0.5. The clus- 
ter order parameter a (thin lines) is expressed in units of ctc, 
the biaxiality order parameter q (thick lines) is expressed in 
units of \/cJd. The temperature function a(r) is scaled by 
the constant a* = 36^/16(7. 
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FIG. 2: Temperature dependence plots of the calculated order 
parameters a (thin lines) and q (thick lines) for fe < 0, A = 
1.05, u = 0.5, for the indicated values of the applied electric 
field E, expressed in units of (a*c^/2e/i^) ^^"^ . 



of the molecules. At constant the transition from 
to Nu'^^ is accompanied by a jump in q that is propor- 
tional to E^. A phase transition of this type, reflected 
on the abrupt change in the biaxial response of a uniax- 
ial nematic phase to an applied electric field, has been 
observed by J-H Lee et al.ld^. 

For weak applied fields, a measure the susceptibility 
of the system to field-induced biaxial ordering is pro- 
vided by the "electro-biaxial" coefficient k 



E=0 



In the Nu phase, the value of this coefficient is fixed to 
fciVu — h/c, while in the Nu'^'^ phase it varies with ct 
(and therefore with temperature) according to fc^cbc) = 

(h/c) (1 — (t/(Tc) ^. Thus, at the transition the electro- 
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biaxial susceptibility undergoes a jump of h/{c{X — 1)). 
This can be quite large in systems for which the A ratio 
is close to one. In this case a weak field could induce 
considerable biaxial order to a uniaxial nematic phase. 

The results regarding the application of an electric field 
should be viewed within the practical limitations that 
the assumed alignment of the director n imposes on the 
possible magnitude of the applied field E: To actually 
maintain the uniform orientation of n along the rubbing 
direction [Z axis) of the plates for a nematic of positive 
dielectric anisotropy, the field strength cannot exceed the 
critical value Ec for the Fredericzs transition that reori- 
ents n along the Y axis. Apparently, no such limitation 
applies for nematics of negative dielectric anisotropy; in 
this respect, such systems would be advantageous for 
the study of electric field-induced biaxial order in the 
iVu'"^'' phase. 

To summarise the results, three nematic phases are 
identified in the absence of an applied field: a purely 
uniaxial phase, a spontaneously biaxial phase and an in- 
termediate, macroscopically uniaxial phase consisting of 
biaxial clusters that are randomly oriented. The appli- 
cation of an electric field could induce substantial biaxial 
order to this intermediate {Nu'^'^) phase. The switching 
of this phase between an optically uniaxial and a biaxial 
state as well as its possible transformation to a nematic 
phase (Nu) in which the application of a weak or mod- 
erate electric field does not induce a measurable optical 
biaxiality, are in agreement with experimental observa- 
tions on bent-core nematics f^. 

The above properties of the Nu'^^ phase are not in con- 
tradiction with the experimental observations of biaxial- 
ity by NMR[3i] and XRD \^ since in both cases an align- 
ing field, magnetic or electric, is present. Furthermore, 
according to the induced nature of biaxiality suggested 
by the present analysis, the measured values of the bi- 
axial order parameter by the two experimental methods 
should in general differ, as in fact they do, because the 
biaxiality-inducing fields and mechanisms are different 
in the two methods. It is also worth noting that the im- 
mergence of macroscopic biaxial nematic ordering from 
the collective alignment (spontaneous or field induced) of 
clusters is supported by the experimental observation of 
biaxial order in nematic tetrapodesfH|. Here, the cova- 
lent lateral grouping of the nematogen components into 
quartets promotes the clustering which, in turn, enhances 
the biaxial tendency relative to that of the non-bonded 
nematogens. 

The free energy expression in eq (|4]) can be readily ex- 
tended to include the possibility of polar ordering within 
the biaxial clusters. In close analogy with the formula- 
tion of the a parameter a parameter, p is introduced to 
describe the average magnitude of polar ordering within 
the clusters in a direction transverse to n. The net trans- 
verse polarity of the sample is quantified by means of a 
vector order parameter pA that couples linearly (paEa) 



to the applied field and quadraticly (paPbQab) to the 
biaxial order parameter q. The additional phases de- 
scribed by the extended Landau- deGennes expansion 
include a macroscopically uniaxial nematic phase of biax- 
ial and polar clusters, N^^'^\ and a polar-biaxial nematic 
phase Npl'"^\ Details on the possible phase transitions 
and field-induced effects are presented in a forthcoming 
publication. 

All the results discussed here are based on a simpli- 
fied formulation of the theory wherein perfect uniaxial 
nematic order is assumed and therefore molecular rota- 
tions are restricted in two dimensions. As a result, any 
dependence on the degree of uniaxial nematic ordering 
is suppressed. Removal of this restriction makes the for- 
mulation more elaborate and modifies the details of the 
phase transitions. However, the essential findings regard- 
ing field-induced biaxiality and polar ordering are pre- 
served. These findings offer new insights into the nature 
of phase biaxiality and the related nematic-ncmatic phase 
transitions and broaden the current views on what could 
be considered as a "biaxial nematic" LC for the pur- 
poses of electro-optic device applications. In particular, 
the possibility that some of the known uniaxial nematics 
could in fact consist of randomly distributed biaxial clus- 
ters, suggests that it might be interesting to study more 
closely the electro-optics of certain "uniaxial" nematics, 
specially those with negative dielectric anisotropy. As the 
possibility of fast, field-induced, switching between uni- 
axial and biaxial (and possibly polar) states provides a 
new concept for the design biaxial nematic devices, such 
studies might also be of practical importance. 
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